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Abstract 

In this paper we study the general scenario of an effective theory coming from the com- 
pactification of a higher dimensional theory in a string inspired setting. This leads to 
gauge coupling unification at an intermediate mass scale. After having computed all the 
threshold corrections (due to Kaluza-Klein modes) to the running of the couplings of the 
MSSM we embark in a detailed phenomenological analysis of the model, based on the nu- 
merical package DarkSUSY, to find constraints on the scenario from Dark Matter data. 
The mass spectrum of the theory does not have tachyons. Moreover we find that the neu- 
tralino is still the LSP with a relic density compatible with the most recent experimental 
data. With respect to the standard mSUGRA scenario we find that the neutralino is hig- 
gsino like in most of the parameter space. Our modifications to the DarkSUSY package 
will be shortly available upon request. 



1 Introduction 



The idea of extra dimensions has brought to an explosion of activity in the field of high 
energy physics phenomenology. The unification of the couplings of the theory at an 
intermediate mass scale is in fact particularly appealing. Narrowing the gap between the 
energy scale that can be probed by current experimental techniques and the unification 
energy scale makes indirect test of the latter more reliable. 

The connection between the supersymmetry breaking energy scale and the size of the 
compactified extra dimensions first appeared in PJ E] . At that time the result was consid- 
ered a negative one: a breaking scale of the order of the TEV implied a compactification 
scale of « 10 18 GeV, a region in which perturbative computations in string theory are 
ruled out. This judgement was influenced by the fact that the fashionable string theory 
at the time was the heterotic string and that the other consistent string theories were 
not considered to be relevant from the phenomenological point of view. It was only later, 
with the discovery of dualities among different string theories [3] that this scenario became 
appealing [3]: all consistent string theories gained in fact equal status. Non-perturbative 
computations for a certain string theory became perturbative ones for its dual. Being in 
the perturbative regime is only one of the desirable features that the theory should have. 
Another major hurdle for a model building directly coming from string theory is given by 
the presence of moduli. These are scalar fields which appear after the compactification 
and give flat directions to the potential of the theory. How to fix this problem is not clear 
yet, even if encouraging progresses have been done recently [SJ 13 El E] • 

Other inputs on the relevance of the extra-dimensions scenario for phenomenology 
could be given by a comparison with experimental data of present and sensitivities of 
future experiments. The approach in which a low energy model is derived directly from a 
string theory is called top-down, while trying to prune models on the basis of low-energy 
data is called bottom-up. Eventually these two approaches will merge, but for the moment 
they are still distinct. 

This paper is an example of the bottom-up approach: as suggested in [TU] a simple 
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way to achieve an earlier unification is to modify the running of the couplings in an 
"exponential" fashion. This can be accomplished by keeping into account the effects of 
extra massive Kaluza-Klein (KK) states. The scenario here is quite generic. A detailed 
knowledge of the compactification mechanism is not needed. At the start only the gauge 
sector and the Higgs fields are allowed to have a tower of KK states. At a later time also 
part of the fermionic sector can have a tower. 

In this paper we analyse the implications of this scenario for the problem of dark 
matter using the numerical package DarkSUSY ^T] together with a modified version of 
the ISASUGRA ^2] package. In order to do this we compute the one loop running of all 
the couplings in the minimal SUSY extension of the standard model (MSSM) in the extra 
dimension scenario of Ref . ^H] . To do this efficiently we use the method of the effective 
potential |14j . See the lecture notes in Ref. for a nice summary. The computation of 
the running of the gauge and Yukawa couplings was performed already in Ref. jTU] , while 
the running of the soft parameters was first studied in JH] using the spurion formalism. 
Instead we used a diagrammatic analysis finding some discrepancies between our results 
and that of the previous reference. 

This is the plan of the paper: in Section |21 we outline the theoretical framework and 
the method of the effective potential. In Section El we present the renormalization group 
equations in presence of extra-dimensions that we implemented in the numerical package 
ISASUGRA. In Section |U we show the results of our computations: the spectrum of the 
theory has no tachyons. The neutralino remains the LSP particle in almost all of the 
parameter space we have explored. If no fermions are allowed to have a KK tower then 
the neutralino is higgsino-like. Otherwise it is bino-like for tan (3 < 30 and higgsino or 
mixed for tan/3 > 30. 
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2 Extra-dimensions, renormalization and the effec- 
tive potential 



In this Section we deal with the theoretical issues which must be faced to adapt the 
DarkSusy package to the extra-dimensions scenario of ^U]. In the first subsection we 
review the main aspects of the extra-dimensions scenario of JU] while in the second 
subsection we review the computation of the one-loop /3-functions of the MSSM with the 
method of the effective potential |14[ I15j. 

2.1 An extra-dimensions scenario 

The MSSM emerges after the extra dimensions are compactified on circles of fixed radius 
R. The energy scale R^ 1 = /x exceeds experimentally attainable energy scales and 
represents the threshold beyond which the effects of extra-dimensions can be felt. Let 5 be 
the number of extra-dimensions. The effect of the compactification of a 4 + 5-dimensional 
field along 5 circles of radius R is to allow for the presence of four- dimensional massive 
excitations, called KK modes, of mass m\ = + n ■ n/ R 2 . n is a ^-dimensional vector 
with rii G Z. mo is the mass of the zero-mode. If an additional Z2 is imposed on the 
compactified dimensions, the circles change into orbifolds. Fields sitting at the orbifold 
fixed points are not allowed to have KK modes. This is a convenient mechanism to 
select which fields have massive KK excitations. Let us now denote by rj the number 
of generations of MSSM chiral fermions. In the simplest possible scenario rj = 0: no 
chiral MSSM fermions have KK excitations. We will discuss scenarios with both 77 = 
and rj ytz 0. Compactifications with S = 1,2 will exhibit N = 2 SUSY. Higher values 
of 5 will result in theories with an even larger number of supersymmetries. The N — 1 
vector multiplet, together with an additional chiral multiplet 1 , are such that their 
KK modes are organized into an iV = 2 multiplet. Finally the two chiral fields, Q, Q, 
containing the two Higgs fields form an iV = 2 hypermultiplet. For 77 > 1 also the chiral 

multiplet containing one of the chiral fermions and its mirror are organized into an N = 2 

lr This chiral multiplet must be odd under the Z2, so that its zero-modes are missing, to recover the 
expected N = 1 MSSM. 
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hypermultiplet. As we said earlier, this N = 2 SUSY for the KK modes is welcome since 
it allows the gauge couplings to be computed a la Wilson. To make the Lagrangian of the 
theory really N = 2 SUSY invariant we must add a term, which couples the chiral 

field in the N = 2 multiplet with the fields in the hypermultiplet. This term is only there 
for the KK modes: zero modes have, in fact, a. N = 1 SUSY. To test how robust is our 
scenario to this term we have performed computations with and without this coupling: as 
we will see when we will discuss our results, there seem to be minor differences between 
the two cases which we will dub Higgs N — 1 and N = 2. 

2.2 Computations of the one loop /3-functions 

For energy scales larger than /xq, KK excitations become possible and their effects have to 
be incorporated in the computation of the /3-functions. Given the large number of such 
functions to be computed in the MSSM we prefer to carry out this computation using the 
method of the effective potential ^3J QSj. In fact, in SUSY theories all divergences are 
found in the wave-function renormalization of chiral multiplets and in the renormaliza- 
tion of gauge coupling constants (equivalent to wave-function renormalization of vector 
multiplets). At the one-loop level, it is then sufficient to consider the renormalization 
of the scalar potential to determine all the renormalization constants of the theory. The 
one-loop divergent contributions to the scalar potential are given by 



where A is a cut-off mass parameter. The first term contains the quadratic divergences. 
StrA4 2 (z) is considered as a function of the scalar fields z a and their vacuum expectation 
values. StrA4 A (z) is the analogous supertrace for the fourth power of the mass matrices. 
The second term contains all logarithmic one-loop divergences. 

The scalar potential is 




(2.1) 



V(z a , z a ) = V susy (z a , z a ) + V soft (z a , z a ) 



(2.2) 



where V susy is the standard supersymmetric potential 



V susy (z a ) z a ) = f a r+ l -D A D 



(2.3) 
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and where we used the standard notation f a = df/dz a , f a = df/dz a . More subscript 
and superscript indices denote higher order derivatives. The D-terms read 

D A = -gz a T Aa b z b (2.4) 

in terms of the generators T A of the gauge group for the representation of the chiral 
multiplets. The superpotential will be written 

/ = \ilabZ a Z h + ]:fa b cZ a Z b Z C (2.5) 

2 o 

Following the notation in ^H], we set z a = H u , Hd, Q, L, u, d, e for the chiral superfields of 
the theory. Since the top quark, bottom quark and tau lepton are the heaviest fermions 
in the Standard Model, we will use the approximation in which only the 33 component of 
the Yukawa matrices is important. Then, in terms of the components with respect to the 
weak SU(2), the relevant fields become Q3 = (t,b),L s = (u t ,t),H u = H®), Hj, = 
(Hj, Hj),Us = t , c?3 = b,es = f and the superpotential reads 

/ d 2 9f a f a « yt(tiH°u ~ ibH+) - y b (tbH^ - bbH° d ) - y T (u T fH d - frH° d ) 

+ rtH+H;-H°H%) (2.6) 

V so f t (z a , z a ) keeps in account soft breaking terms. These are gauge invariant terms which 
do not add quadratic divergences that is StrM. 2 (z) does not receive any new field depen- 
dent contributions from these terms. Keeping in account gauge invariance we only need 
to add scalar and fermionic terms 

V soft (z a , = (m 2 ) a b z a z b + + £(z) (2.7) 

where £ is an arbitrary gauge invariant cubic function of z a 

£(z) = i [1)a z a + i {2)ab z a z b + t mahc z a z b z c (2.8) 

V so f t (z a , z a ) only adds constant terms to the squared mass matrix and ()2.7|) is the most 
general expression for scalar soft breaking terms. They contain scalar mass terms of the 
form zz, z 2 + z 2 and trilinear analytic interactions z 3 + z 3 . The spin 1/2 mass matrix 
instead 

——A AB \ A X B - -A AB \ A X B (2.9) 
2 2 



will only give a constant contribution A AB A BA to trMy 2 . This term is then a soft term 
too. All other additions to M\j 2 are not soft in general. These soft breaking terms contain 
all we need to find an acceptable spectrum of SUSY particles. We have soft terms to raise 
the masses of scalar quarks and leptons, and we also have gaugino mass terms to raise 
the mass of gluinos and photino. Using the same notation as above, denoting by a t , a b , a T 
the 33 component of the matrices £(3) u qh u , £,(3)dQH d , £,(3)eLH d and assuming that the squark 
and lepton matrices are flavour blind we get 

1 



V f 



soft 



2 {A 66 gg + A 22 WW + ABB) + c.c. 



m 2 Q J*t - m 2 Q3 b*b 



- m L v T v T - m L T t — m u t t — m d b b — m g T r 

- a t {ttH° u - ibH+) - a b {tbH d - bbH° d ) - a T (u T fH d - frH Q d ) 
+ [bWH- - H° u H° d ) + c.c] - m\ u [(H:)*H: + (H° u yH°] 

- ™ 2 H d [{H~)*H- + (H°)*H°] (2.10) 

g, W, B are the gaugino, wino and bino. Let us now compute the supertrace of the quartic 
mass matrix for arbitrary values of the fields z a . Let's start with the mass matrix for spin 
one 

The fermionic mass matrix has instead the form 

Ml/2 ={-iV2D*« A AB ) (2 ' 12) 
From this matrix one obtains 

-2trM\ /2 = -2f ab f bc f cd f da -lQf ab f ac D A D Ab - lQ(D A D Ab )(D Ba D B ) 
+8f ab D Aa A AB D Bb + f ab D A A AB D b - \<6A AC A CB D A D Ba 

_ 2A AB A BC A CD A DA 

The mass matrix squared for the scalar fields is given, as usual, in terms of second deriva- 
tives of the potential 

' V b V ac 



where 

V ab = fabcf c + D A D A + Vab (2.15) 

yab = fobcj c + D Aa D Ab + -ab ^ ^ 

Vb = f ac fbc + D Aa b D A + D Aa Dt + {m 2 )l (2.17) 

Note that the mass matrix squared is field independent. The soft terms give 

trM 2 soft = 2(m 2 y a - 2A AB A BA (2.18) 

which is also field independent. 

For the trace of the fourth power of the mass matrix, after some calculations we finally 
obtain 

StrM A = +2g 2 [T{R) + 2C{R) - 3C(G)]D A D A - 8g 2 C(R)f a f a + 

+4f ac f bc D Ab a D A + 2f abc f abd f d f c - 2g 2 C(R)[z a ff ab + z a f b f ab ] + 
-8g 2 C(R)A A {f a z a + f a z a ) - 16g 2 C(R)A 2 A z a z a - 2A\ + 
+2f abc f cVab + 2f abc f c V ab - 2g 2 C(R)(r] a z a + rj a z a ) + 2«° 6 + 
+4f ab f ac (m 2 y b + 4g 2 C(R)(m 2 y b z a z b + 2(m 2 ) a b (m 2 ) b a + 
+A(m 2 ) a b (D A ) b a D A (2.19) 

From (J2.19)) we can compute the coefficients of the terms with different powers of the 
variables z, z. The renormalization of the gauge coupling constant g and the Yukawa 
couplings f abc is not affected by the soft terms which do not contribute to the dimension 
four terms zzzz. Then one defines the following renormalization constants 

5 a b - U) z b (2.20) 



with 



9 = 0-+p)9 (2-21) 

1 ; 1 / 1 / / / / 

fabc = fabc + -t-a fa'bc + ~^b fab'c + fabc' + fabc ° fa'b'c' (2.22) 

flab = P-ab + ^ e a Pa'b + ^b Hab' + Hab Va'b' (2.23) 



p = -2kg 2 [T(R) - 3C(G)} (2.24) 
e a b = Akg 2 C(R)5 a b - 2kX^ (2.25) 



where X 6 a = f acd f bc d and 



Furthermore 



^ = 2kg*[C(A) + C(B)}6 a :6« (2.26) 
f£ d = 2kg 2 [C{A)+C{B) + C{CW:5t5i (2.27) 



fc = ^ln(A/Q) (2.28) 



where A is an ultraviolet cutoff and Q is the energy scale at which the couplings are 
computed. As it is well-known in SUSY theories the couplings renormalize as 

^ = {Z^r^t^y (2.29) 
f abc = {Z^)i{Z^)l{Z^)if a , blc , (2.30) 



From ()2.30|) we infer 



[Z 1/2 ) a b = 5 a b + Akg 2 C{A)5 a b - kX£ (2.31) 



It is now easy to compute the renormalizations of the couplings and masses of the theory 
and from this their runnings. In Appendix A we collect, for the reader's convenience, all 
these results following [15] . It is now easy to specify these formulae to the fields of the 
MSSM and recover the results of ^B] for the running of the MSSM couplings and masses. 

In presence of KK modes the effect of the thresholds is to contribute to the running 
with power-like terms JU] • The general form of the running of a coupling g will be of the 
type 

Q% = {b-b)hi^ + X s b(-) (2-32) 

Xs is a constant depending on the number of compactified dimensions, b is the standard 
one loop value of the coefficient, while b is the contribution of the thresholds. 



3 Running of the couplings 

In this Section we list the running of the couplings of the theory. To find them we used 
the formulae of Appendix A from which we extracted the contributions of the various 
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KK towers. In certain cases we also checked these results against the appropriate Feyn- 
man diagrams which renormalize the theory. Appendix B contains some details of the 
computation which will be fully discussed in [T7j . 

The running of the gauge and Yukawa couplings were first studied in in this case 
the coefficients b a take into account the Higgs and gauge bosons KK tower and, separately, 
the contribution of rj families of chiral fermions. The one-loop coefficients (|2.32|) for the 
running of the gauge couplings are 

33 

b a = (-,1,-3) 
5 

b a = (3/5, -3, -6) +77(4, 4, 4) (3.1) 

where a = 1,2,3 labels the three gauge groups £7(1), SU(2), SU(3). The £7(1) generator 
is rescaled by a factor 3/5. 

As we said earlier, it is convenient to assemble the two N = 1 Higgs chiral fields 
into an N = 2 hypermultiplet. To be invariant under N = 2 SUSY, besides the kinetic 
terms for the vector multiplet and the hypermultiplet the lagrangian must also contain 
an interaction term of the following form 

J d 2 6 (V2H u $H d + fiH u H d ^ + h.c. (3.2) 

where /1 plays the role of the usual /i-term appearing in the MSSM while <£> is the chiral 
multiplet contained in the N = 2 vector multiplet. These are the only admissible terms 
for an N = 2 theory. The first term is a coupling between the hypermultiplet and the 
vector multiplet, and it is needed to enlarge the amount of supersymmetry. The one- loop 
diagram built from the interaction vertex associated to this term exactly cancels the one 
loop diagram of Fig. 6d of JU] ( m whose loop KK states can propagate) thus ensuring 
that the Higgs wavefunction is immune to the effects of the extra space time dimensions. 
The first term in ()3.2)1 is not present at the level of zero modes because $ is odd under the 
action of the Z2 of the orbifold. In this way even if the Higgs zero-modes are assembled 
into an iV = 2 hypermultiplet they still have an anomalous dimension. 

To check how robust the scenario is with respect to the presence of this term, we have 
studied two different cases which we call Higgs N = 1 and N = 2. In the former case we 
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assume that the KK states behave as their zero-modes and that the first term in ()3.2|) is 
absent. In the latter case this term is switched on. 



3.1 Higgs N = 1 

3.1.1 Couplings in the SUSY Lagrangian 

In this case the anomalous dimension receives a contribution from the diagrams in which 
particles from the KK tower of the gauge bosons circulate in the loop leading to a power- 
law contribution. The wavefunction renormalization factor of the field % is 



Zi = l In 7-7. (/4o) 

Z7T yU 



(-) -i 



leading, for the Higgs fields to 



(3.3) 



1h = cti ao 

' " 10 2 

3 3 

7ff. = ai 0:2 

' d 10 2 

with ctj = gf/4:iT. For what the Yukawa couplings are concerned, we get 

6 



2 G? 
167f df * 



2/t 



3|^| 2 + X^^ (3|^| 2 + |y 6 | 



16 



"53 



9 1 3 o 

-3<7 2 - - 9l 



at 



Ub 



3\y b \ 2 +\y T \ 2 + X 5 (^j (j^| 2 + 3|y b | 



1ft 2 rf 
16n dt y 



16 2 2 ^ 2 

-y£ 3 -3^2 - Y5^i 



\Vt\ 2 + 3\y b \ 2 



V/W 

Finally we analyze the (supersymmetric) Higgs mass 

d 



1671 IS 



(3||/ t | 2 + 3|i/ fe | 2 + |i/ r | 2 ) 



o 2 3 2 

-3# 2 - 



(3.4) 
(3.5) 



(3.6) 



(3.7) 



(3.8) 



(3.9) 
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In the Higgs N = 2 case fi will not be renormalized as a consequence of a non-renormalization 
theorem. 



3.1.2 Soft terms in the minimal scenario 



While the effective potential method is powerful in determining the runnings for the 
parameters in the SUSY Lagrangian, it is difficult to apply it to the computation of the 
soft terms in presence of the KK states. This is why we have decided to also carry on 
an analysis in terms of Feynman diagrams. In Appendix A we then collect the general 
form of the running of the couplings and masses appearing in (J2.2J) . In Appendix B we 
collect all the contributions to the renormalization of the couplings and masses in ()2.10|) . 
As a check, the interested reader can sum (jB.lj) - ()B.4|) to recover the result of (jA.5|) when 
specialized to the parameter b in ()2.10|) . The computation of the effect of extra dimensions 
is now straightforward: each time that in a Feynman diagram KK states are allowed to 
circulate in the loop (external states can only be zero modes), the contribution of the 
diagram must be added to the b coefficient in 1)2.32)1 . We finally remind the reader that 
in an orbifold compactification the KK momentum is conserved only in the bulk. The 
wavefunction in the compactified dimension must, in fact, be expanded in a basis which 
is invariant under the discrete group acting on the compactified dimension. Furthermore 
the orbifold fixed points break translation invariance along the extra dimensions. 

Following the previous prescription we obtain the /5-functions for the trilinear couplings 



2 d 
167r — a t 

at 



2 d 
I6n —a b 

at 



9a t \y t \ 2 + X s ( 
+2a b y* b y t + y t 



32 



a* 



16 

%t| 2 + N 2 -yS 3 2 -32 2 2 



26 



13 
15 



9i 



3 g 2 3 M 3 + Qg 2 2 M 2 + — 



9a b \y b \ 2 + a b \y T \ 2 + 2a T y*y b + X s ( 



A 
Q~ 



a b 



9M 2 + M 2 + 



-y#3 - 3# 2 - —g\ + 2a t y*y b + 



+Vb 



32 14 

— g 2 3 M 3 + 6g 2 2 M 2 + -^g\M x 



(3.10) 



(3.11) 
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16n 2 ^a T = 3a T \y b \ 2 + 3a T \y T \ 2 + X s (-^-) S { a T 9\y T \ 2 + 



o ^9 



+ Ga h yly T + y 7 



6g 2 2 M 2 + —g\M x 
o 



(3.12) 



The /9-function for the b parameter (crucial for the electroweak symmetry breaking mech- 
anism), is given by 



lQn 2 ^-b = b(3\y t \ 2 + 3\y b \ 2 +\y T \ 2 )+fi[Qa ty ; + Qa b y* b +2a T y: 



dt 



(3.13) 



+X 5 (J^j lb {-3g 2 2 - ^ + fi (+6g 2 2 M 2 + \g\M^j 



Now we can show the /3-function for all the soft squared masses. For the two Higgs mass 
parameters we have 



\^ 2 jm 2 Hu = Q\y t \ 2 (m 2 Hu +m 2 t +ml)+Q\a t \ 2 + ^g 2 tr'(Ym 2 ) 

+Xs (A)' {-Gg 2 2 \M 2 \ 2 - \gl\M.\ 2 + \g\ « - m 2 Hd ) 



lo7r —m M 
dt Hd 



%b| 2 (m 2 Hd + m 2 + mf) + 6|a fe | 2 + 2\a T \ 2 
+2\y T \ 2 (m 2 Hd +m 2 L + w 2 ) - -g 2 tr\Ym 2 ) 

+Xs (A)' {-Q 9 2 2 \M 2 \ 2 - \gl\ Ml \ 2 - \g\ (m 2 Hu - m\ d ) 



(3.14) 



(3.15) 



while the equations for the squarks mass terms read 
2\y t \ 2 m\ + 2\y b \ 2 m\ + \g\tr\Ym 2 ) 



i r 2 ^ 2 
l07T — m f 

dt * 



(3.16) 



+Xs (A)' | - f S 3 2 |M 3 | 2 - 6^|M 2 | 2 " ^imi 2 + 2\a t \ 2 + 2\a b [ 



1 r 2 ^ 2 
l07T — m f 



+2\y t \ 2 (m 2 Hu + m 2 ) + 2\y b \ 2 (m 2 Hd + m 2 ) + ^g\ (m 2 Hu - m 2 Hd ) 
A\y t \ 2 m 2 t - yitrXYm 2 ) 

+X S (A)' | _ ^|M 3 | 2 - ^IM,! 2 + 4h| 2 
+4|y t | 2 (m 2 ^ + m 2 ) - -g\ (m 2 Hu - m\ d ) 



(3.17) 
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167r 2 |m? = ^ytfmt + pltr'iYm 2 ) (3.18) 

+ ^(^)^f^ |M3|2 -^ |Ml|2 + 4N2 
+4M 2 (m 2 Hd + ml) + -g\ {m\ u - m 2 Hd ) 

Finally, for the sleptons masses we have 

IQ^^-ml = 2\y T \ 2 ml-\g[tr\Ym 2 ) (3.19) 
at 5 

+X S | _ 6^ 2 |M 2 | 2 - ^IM^ 2 + 2|a r | 2 

+2|y T | 2 (m 2 ^ + m|) - ^gl (m 2 Hu - m^J 

lGTr^ml = ^\y T \ 2 m 2 L ^\g\tr\Ym 2 ) (3.20) 
at 5 

+ - Y ^) f {-ir J?|Ml|2+4K|2 

+% T | 2 (m^ + ml) + j^ 2 (m 2 ^ - m 2 Hd ) 

where, in order to avoid cumbersome expressions, we have introduced the quantity 

tr\Ym 2 ) = tr(Ym 2 ) - m 2 Hu + m\ d (3.21) 

where the usual trace over the hypercharge, defined in the MSSM, is given by 

tr(Ym 2 ) = m 2 Hu -m 2 Hd +m 2 Q -2{ml + ml + mf) (3.22) 
+m 2 j + m\ + ml - m 2 L + ml + m| + m 2 

In our computations we have assumed the first two families to be degenerate in mass, 
such that 



m Q 


= 2rn 2 Ql 


+ m Q 3 


m 2 L 


= 2m 2 Li 


+ ml i 


ml 


= m l 






= m 2 5 





(3.23) 
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3.1.3 The non minimal scenario 



In the non minimal scenario we introduce KK towers also for the matter fields. We must 
distinguish among three different cases: 77 = 1 in which only the third family have a KK 
tower, r) = 2 in which only the first two families have a KK tower and rj = 3 in which all 
the families have KK tower. The choice for the 77 = 2 case is dictated by the usual third 
family approximation and by a constraint imposed by the ISASUGRA code concerning 
the mass degeneracy of the first two families. 

One of the consequences of having KK towers for the matter fields is that they no 
longer live at the orbifold fixed points and they are allowed to be in the bulk. 

The anomalous dimensions of the Higgs fields are the same in the 77 = and 77 = 2 
case. For 77 = 1, 3 we find 

1h u = 1h u (3.24) 
lH d = lu d (3.25) 

For the chiral fields t, b, r we have a similar situation and the coefficients of the /5-function 
have only a power-law contribution. 

Let us now discuss the soft terms: for the trilinear couplings and the b parameter, we 
must distinguish between the 77 = 2 and 77 = 1, 3 cases. In the former case the running is 
the same of the minimal case due to the third family approximation. In the latter case, 
once again, the /3-functions coefficients have only a power law behaviour. 

Finally the running of soft scalar masses in the non minimal cases is given by a 
logarithmic term multiplied by 

^(Ym 2 ) (3.26) 

and a power law term given by 

tr(Ym 2 ) - tr(Ym 2 ) (3.27) 
In the 77 = 1 case the primed trace is defined as 

tr' v=1 (Ym 2 ) = 2m 2 Qi - 2 (m? + mf) + m\ + m\ - 2m 2 Li + 2m? (3.28) 
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while for r] = 3 



tr' v=3 (Ym 2 )=0 



(3.29) 



because all the families have KK tower and there is no logarithmic contribution. Finally 
for the case rj = 2 we have 



^ =2 (Fm 2 ) = m 2 Q3 - 2 (mf) + mf - m| 3 + m? 



(3.30) 



3.2 Higgs N=2 

With respect to the previous Section, in this case the changes are given by the presence 
of the term ()3.2j) . Following ^0] the wave function renormalization of the Higgs fields are 
given by 



lH d = 



(3.31) 
(3.32) 



The one-loop beta functions for the Yukawa couplings are then 

uU \\ ( — ) 

10 ai 2 



!6vr —y t = y t 



dt 



3M 2 - Tn9\ ~ hi + X S (A) (3|^| 2 + M 2 - ^g! 



3 2 _ 17 2 
2 #2 3Q 9i 



(3.33) 



2 ^ 

16vr — y 6 = y 6 



16 2 3 2 1 2 



(3.34) 



IK 2 d 

16vr — y r = y T 
dt 



^2 ^ 2 

2^ 2 ~ 2 dl 



(3.35) 



For the (supersymmetric) Higgs mass, in virtue of the non-renormalization theorem, we 
have 

1" = < 3 - 36 > 
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This equation holds for any 77. The /^-functions for the b parameter reads 
d. 



16tt z 



dt 



b (3\y t \ 2 + 3\y b \ 2 + \y T \ 2 ) + fi (Qa t y* t + Qa b y* h + 2a T y* T ) 



(3.37) 



3.2.1 Soft terms in the minimal scenario 



We begin with the trilinear terms, which look like 



167T — at = 
at 



if? 2 d 
167r — a b 

dt 



a«(%«l 2 -^-^) + x 5 (A )Na< 

+a t (9\y t \ 2 + \y b \ 2 ) + 2a b y* b y t 
+Vt 



16 2 _ 3 2 _ 17 2 
• 3 #3 2 #2 3Q ^i 



32 26 
y# 2 M 3 + 6(7 2 2 M 2 + y^iMx 



a& (% 6 | 2 + \y T \ 2 - ^gl - + 2a T y*y b 



a b 



1 ,9 1 ,9 16 9 ^2 "^2 
%&l + |yt| - y# 3 - 2^2 - g0i 



+2/6 



32 , 9 14 , ' 

y ^M 3 + 6^M 2 + ^1^1 



167r 2 ^a T = a T ( 3 1 1 2 + 3\y T \ 2 - ^gf - ) + (wi,!)l!) T 



%r| 2 - - 7^1 



(3.38) 



(3.39) 



6^2 2 M 2 + ^s?Mi 




(3.40) 



The /3-functions for the soft squared masses of the Higgs fields are 



1 r 2 d 2 
107T -m H 



d 

— T lib 1 1 

dt Hd 



16n 2 ^m 2 



= 6|t/ t | 2 (m% u + m 2 + ml) + 6\a t \ 2 + -gltr{Ym 2 ) (3.41) 

-6<? 2 2 |M 2 | 2 -\gl\M.\ 2 
5 

= Q\y b \ 2 {m 2 Hd +m 2 + ml)+6\a b \ 2 + 2\a T \ 2 (3.42) 
+2|y T | 2 (m 2 Hd + m\ + m%) - -gltr{Ym 2 ) 
-6( ? 2 |M 2 | 2 -^ 1 2 |M 1 | 2 
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while for the squarks we have 

167r 2 ^m t 2 = 2\y t \ 2 m 2 l + 2\y b \ 2 ml+ l -glti J {Ym 2 ) (3.43) 

A V f 32 



A', \q-J \~ ySsl^r - 6^|M 2 | 2 - -gilM^ 2 + 2\a t \ 2 + 2\a b \ 

+2\y t \ 2 (m 2 Hu + m 2 ) + 2\y b \ 2 (m 2 Hd + m 2 ) + X -g\ (m 2 Hu - m 2 Hd ) 

d 4 
16tt 2 — m\ = 4\y t \ 2 m 2 - -gH'r'iYm 2 ) (3.44) 
at 5 



+4|y t | 2 (m^ + m 2 ) - ~# 2 (m^ - m|J 



d 2 
16tt 2 — m 2 = 4|y b | 2 m 2 + -gHr 1 (Ym 2 ) (3.45) 
at 5 

+A 5 | - f 9lW,\ 2 - ^|Mx| 2 + 4|a,| 2 

2 

+4M 2 (m^ d + m 2 ) + -g\ (m 2 Hu - m|J 
Finally the equations for the leptons read 

16tt 2 — mf = 2\y T \ 2 m 2 f -^gftr(Ym 2 ) (3.46) 

+X S (A)' | _ 6s 2 |M 2 | 2 - ^ 2 |M X | 2 + 2|a T | 2 

+2|y T | 2 (m 2 ^ + - jj(7 2 (m 2 ^ - m^J 
167r 2 A m 2 = 4\y T ^ m l + ^ g l tr '(Ym 2 ) (3.47) 

+X '(^) i {-^I M 'I 2+4 ^I 2 

+4|y r | 2 (m 2 ^ + m 2 ) + j^ 2 (m 2 ^ - m|J 

Our results (|3.43)) - (j3.47)) for the masses of the squarks and sleptons differ from those 
reported in jTH] in two respects: on the one hand the terms proportional to Y and m 2 H — 
m 2 H d i due to the U (1) gauge factor 2 , are absent in |XH|. On the other hand in (j3.43|) - (j3.47|) 



2 For an explicit evaluation of this term, given by the diagram (|ij.l4fl . see for example (4.12) in ^j] 
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there are terms coming from the diagrams (|B.13|) which do not get contributions from the 
KK states and therefore do not have a power-law type running. 



3.2.2 Non minimal scenario 

The results we found in Section l3.1.3l for N = 1 are still valid. The only difference concerns 
the soft mass terms: for the Higgs masses, rn 2 Hu , no diagrams will contribute to the 
power law. The equations are then the same as in the MSSM. On the other hand the 
squark masses have only power law contributions. For the other soft terms the equations 
are the same of the minimal scenario, regardless of the value of n. 

4 Phenomenology 

In this Section we present the results for the low energy phenomenology of our scenario. 

To compute the weak scale parameters we used the publicly available codes ISAS- 
UGRA [HI and DarkSUSY [1 1J . We modified the ISA JET routines in order to take into 
account the power-law running of the gauge couplings, of the Yukawa couplings and of all 
the soft terms following the procedure outlined in the previous Sections. We solved the 
differential equations given by the running of the couplings and masses imposing bound- 
ary conditions as in the usual mSUGRA scenario ^3]- In this way we have, at the new 
unification scale M' GUT , a universal scalar mass mo, a universal gaugino mass and a 
common trilinear coupling Aq. The only other parameters which need to be specified are 
tan (3 and sign(fi). 

For a fixed choice of the extra-dimensional parameters /io, 5, rj we performed a detailed 
scan of the (mo,mi/ 2 ) parameter space keeping fixed tan/3, A and sign(fi). We computed 
for every model the thermal relic density 




(4.1) 



dn 



'X 



+ 3Hn 



(a A v) [{nf - [n eq f] 



(4.2) 



dt 



'X 
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where H = a/a is the Hubble expansion parameter (with dot denoting time derivative), 
a = a(t) is the scale factor of the Universe, (<Jav) is the thermal averaged annihilation 
cross section and n eq is the equilibrium number density. We considered only tree-level 
annihilation processes in the computation of (<jav). This implies that we do not have 
to worry about the presence of the KK states. In fact KK modes are too heavy to be 
produced in the intermediate tree-level annihilation channels. This is guaranteed by the 
fact that the neutralino mass (i.e. the mass of the lightest supersymmetric particle) is at 
most of order of 1 TeV for every model we considered. For the relic density computation 
we have taken into account all the possible coannihilations with other sparticles. Because 
the mass of the first KK excited level is of order p there is no possible coannihilation 
between the neutralino (a 0-mode particle) and the KK particles. In this way the standard 
treatment of the coannihilations is still valid. The possible difference with the standard 
scenario may arise from the Hubble expansion rate H . In the standard framework, this 
term is responsible for the density n x oc a -3 , in absence of interactions. Unfortunately 
in the literature there is not a general expression of the Hubble expansion rate H for an 
arbitrary number of extra dimensions 5. 3 There are some interesting issues concerning the 
D = 6 dimensional case in [27JI2E] and also in [231110]. We considered explicitly only the 
case D = 5 |3T] and we argued that the conclusion still holds in the higher dimensional 
case jS2]- Let us consider the case 5 = 1. Following [20] and [33] we have 



with M 5 being the five dimensional Planck mass. The third term, where C is an integration 

constant, is referred to as the dark radiation term. The second term proportional to p 2 

and the dark radiation term are the new ingredients in the brane world cosmology and 

they lead to a non standard expansion law. The second term is the dominant one when 

p 3> Pq. In order to satisfy this condition we must have M 5 < 10 4 TeV. 
3 See O ESI EH] for general reviews on the subject. 




(4.3) 



where 



p = 96ttGM| 



(4.4) 
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The D = 5 Planck mass is related to the D = 4 Planck mass through the relation 

(M 5 ) 3 R = M 2 P (4.5) 

where R = l/po is the compactification scale. In our model we obtain the values collected 
in Table ^ This result holds in general in model with flat extra space-time dimensions. 
Typical values of the halo density p are of order of 1 GeV/cm 3 so the p 2 term is always 
subdominant. Hence the usual Friedmann equation holds. 



(M) = l/R 


m 5 = (p M 2 P y/i 


po = 967rGM 5 6 


(Gev) 


(GeV) 


(GeV/cm 3 ) 


10 5 


2.15 • 10 14 


10 8 


10 6 


4.64 • 10 14 


10 10 


10 8 


2.15 • 10 15 


10 14 


10 10 


1 • 10 16 


10 18 


10 12 


4.64 ■ 10 16 


10 22 



Table V. Values of the D = 5 Planck mass M 5 and of po for different values of the 
compactification scale p . 



Let us fix the notation defining the lightest neutralino as the linear combination 

X? = N n B + N 12 W + N n H u + N u H d (4.6) 

where B and W are the bino and wino fields while H u and are the two higgsinos. We 
also define the gaugino fraction as 

Z s = |i\r u | 2 + |i\r 12 | 2 (4.7) 

We say that a neutralino is gaugino-like (in particular in our case bino-like) if Z g > 0.9 
while is higgsino-like when Z g < 0.1. In all the intermediate cases we denote the neutralino 
as mixed-like. 

As explained in SecE] we considered two different scenarios: one in which the two 
zero-mode Higgs fields are chiral N = 1 superfields as in the MSSM and the other in 
which they form an iV = 2 matter hypermultiplet. 
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We first analyze the N = 1 Higgs case in the minimal scenario 77 = in which the 
matter fermions do not have KK towers. We performed a detailed scan in the parameter 
space (m , 7711/2) fixing all the other parameters. In Fig. Qwe show the regions already 
excluded due to either theoretical or experimental reasons for a given model (fio = 10 5 
GeV, 5=1 and r\ = 0). The dark gray region is excluded because the f is the LSP 
rather than the neutralino, the red region is excluded because the models do not achieve 
electroweak symmetry breaking (EWSB) while the blue region is excluded because the 
models do not satisfy the current accelerator bounds (limits on the chargino masses, 
b — > 57, etc.). In almost all the parameter space the neutralino is still the LSP. One of 

fi = W 5 GeV, 6=1, ri=Q, tan/3=10, sgn(/i)=+l, A =0 




500 1000 1500 2000 2500 
mo (GeV) 



Figure 1: Excluded regions in the plane (m , rrii/ 2 ) 

the main result is that, unlike the standard mSUGRA case, the neutralino is no longer 
bino-like but it tends to be a very pure higgsino (see Fig. |2J). This conclusion strongly 
depends from the value of the compactification scale /io- In fact the regions in which the 
neutralino is higgsino-like get smaller for higher values of /io (see the left panel of Fig. El 
for the case /_i = 10 5 GeV and Fig. |U for the two cases /i = 10 8 GeV and ji Q = 10 10 
GeV) . For higher values of /io the higgsino-like region approaches the region excluded due 
to an incorrect EWSB. The previous result is only slightly dependent from the number 
of extra-dimensions 5 as can be seen by comparing the two panels of Fig. |21 The only 
difference is a change in the shape of the excluded regions. In other words the crucial 
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/j =10 GeV, 5=1,7;=0, tan/3=10, sgn(/i)=+l, A =0 £<o=10 GeV, 6=4, i;=0, tan/3=10, sgn(/j)=+l, A =0 




500 1000 1500 2000 2500 500 1000 1500 2000 2500 

m (GeV) m (GeV) 



Figure 2: Contour plots of the gaugino fraction in the plane (m , m^) for different values 
of 5. 

property of the /3-functions is the power-law behaviour rather than the effective power-law 
index, i.e. 5. 

In every contour plot for the gaugino fraction we have shown the cosmologically allowed 
regions. The red regions are those for which the relic density flh 2 satisfies the WMAP 
constraints jHHES] 

0.09 < Qh 2 < 0.13 (4.8) 
while the green regions denote the pre- WMAP constraints 

0.13 < Vth 2 < 0.30 (4.9) 

The cosmologically allowed regions have a huge overlap with the pure higgsino-like region, 
especially for low values of /xo, i.e. /Jo < 10 8 GeV. For higher values of /Jo the cosmologically 
allowed regions tend in general to overlap with a mixed-like neutralino region. 

We also present in Fig. |3] the isomass contour plots for the neutralino for increasing 
/Jo and for fixed 5. As in the case of the gaugino fraction we show the cosmologically 
allowed regions. For higher values of /Jo high neutralino mass contours (around 1 TeV) 
are shifted upwards so that models that possess the right relic density have m x < 500 GeV. 
Increasing the value of Aq, for example A$ = 2500, implies a shift of the cosmologically 
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allowed regions towards higher values of the neutralino mass of about m x > 800 GeV and 
a wider region excluded by the accelerator bounds. Moreover, in this case, the neutralino 
is a very pure higghsino Z g < 0.1 in all the allowed parameter space. 



^o = 10 GeV, (5=1,^=0, tanS=10,sgn(/j)=+l,A =0 /i = 10 s GeV, tf=l, i]=0, tan/3=10, sgn(//)=+l, A o =0 




500 1000 1500 2000 2500 3000 500 1000 1500 2000 2500 3000 

mo (GeV) mo (GeV) 



Figure 3: Contour plots of the neutralino mass (in GeV) in the plane (m , mi/ 2 ) for 
increasing values of /io- 

For increasing values of tan j3 we obtain the same shift of the cosmologically allowed 
regions as in the case of high Aq. Moreover the excluded region in which the neutralino 
is not the LSP grows towards higher values of mo- 

Let us now analyze the non minimal models with rj > 0. As explained in subsec- 
tion we considered three different cases: i] = 1 in which only the third family have a 
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500 1000 1500 2000 2500 3000 500 1000 1500 2000 2500 

mo (GeV) m (GeV) 



Figure 4: Contour plots of the gaugino fraction in the plane (m , 7711/2) for increasing 
values of /io- 

KK tower, rj = 2 in which only the first two families have a KK tower and rj = 3 in which 
all the families have KK tower. 

For rj = 2 we must have /io > 10 8 GeV in order to avoid non perturbative gauge 
couplings at the unification scale (see for example Fig. 4 of JU|). I n this case regions with 
m o ^ 3000 GeV are excluded due to an incorrect EWSB while for tan/? < 30 there is 
no region in which the neutralino is not the LSP. For higher tan/5 the parameter space 
develop a region, for small mo, excluded by EWSB constraints and by the stau being the 
LSP. In almost all the allowed region the neutralino turns out to be a very pure bino (as in 
mSUGRA) though there is still a significant overlap between the cosmologically favoured 
region and the region in which the neutralino is higgsino or mixed-like. The possible values 
of the neutralino mass are lower than in the minimal scenario, with typically m x < 500 
GeV. 

In the case 77 = 1 in which only the third family has KK tower almost all the parameter 
space is allowed for low tan (3. There are only two small regions excluded by the accelerator 
bounds and by the EWSB constraints. The neutralino turns out to be always a bino. The 
cosmologically allowed region is very small and in correspondence with a low neutralino 
mass region, i.e. m x < 50 GeV. For higher tan (3 we have the same behaviour except for 
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the presence of the two excluded regions (EWSB constraints and the stau being the LSP) 
for small ttiq. 




500 1000 1500 2000 2500 3000 2 °° 400 600 8 °° 1200 1400 1600 

m (GeV) "'o (GeV) 



Figure 5: Gaugino fraction in the plane (mo, mi^)- Left panel: non minimal scenario 
i] = 1. Right panel: non minimal scenario rj = 2. 

Finally we consider the r\ = 3 case in which all the families possess KK tower. In these 
models the region for which < 750 GeV is always excluded for any value of tan/3. 
For tan /3 < 30 there is only a very small region in which the neutralino is not the LSP, 
namely for m — GeV and — GeV. For higher values of tan/3 there is another 
region in which the neutralino is not the LSP together with a region excluded by the 
EWSB on the left of the parameter space. The neutralino mass for this kind of models 
is lower, m x < 200 GeV , with respect to the minimal case (left panel of Fig. EJ) . The 
neutralino composition is essentially that of a pure bino except for a region close to the 
excluded regions in which Z g < 0.1 (right panel of Fig. EJ). Once again the relic density is 
in the right range in the regions where the neutralino is higgsino or mixed-like. 

In general non minimal models are disfavoured from the point of view of the cosmo- 
logical relic density because the allowed regions are very small. 

In the scenario in which the Higgs form an N = 2 matter hypermultiplet the results 
about the neutralino composition and the corresponding relic density remain essentially 
unchanged. However there are some differences in the shape of the allowed and excluded 
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/i = 10 u GeV, 5=1,17=3, tan/?=50, sgn(//)=+l, Ao=0 jUo = 10" GeV, *=1, >7=3, tan8=50, sga(jx)=+l, A =0 



mi/2 
(GeV) 




500 1000 1500 2000 2500 3000 500 1000 1500 2000 2500 3000 

mo (GeV) mo (GeV) 



Figure 6: Left panel: neutralino isomass contours (m , m^) for the non minimal scenario 
77 = 3. Right panel: gaugino fraction in the plane (mo, JTO1/2) for the non minimal scenario 
77 = 3. 

regions in the parameter space. For example in the minimal scenario 77 = (see Fig. [7(a)) 
the regions excluded due to EWSB and to the accelerator bounds are smaller than in 
N = 1 case. This holds for all tan/3 and A . Keeping fixed all the other parameters there 
are in general allowed regions in the right part of the parameter space with mo > 2500 
GeV. The cosmologically allowed regions get shifted in order to "follow" the regions in 
which Z g < 0.1. In the non minimal cases 77 > the behaviour is the opposite and the 
parameter space is in general more constrained, although there are cases in which the 
region excluded by the neutralino not being the LSP is absent or at least very small (see 
Fig.|7ic),(d)). 

5 Conclusions 

In this paper we presented a general scenario that involves the presence of large extra 
space-time dimensions. We computed all the threshold corrections (due to KK modes) to 
the running of the couplings of the MSSM, including all the soft terms. Then we performed 
a detailed phenomenological analysis of the model. We found that the neutralino is 
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still the LSP with a relic density compatible with the most recent WMAP data. The 
main difference between our model and the well studied mSUGRA scenario is that the 
neutralino is higgsino like in most of the parameter space. Moreover the mass spectrum 
of the theory does not contain tachyons. These results remain essentially the same in the 
two cases we considered, N — 1 and N = 2 Higgs. 
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(a) Excluded regions in the plane (mo, 
mi/ 2 ) for the minimal scenario r\ = 0. 



/j = 10 5 GeV, (5=1, i]=l, tan/?=10, sgn(/z)=+l, A =0 
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(b) Gaugino fraction in the plane (mo, 
m 1 / 2 ) for the minimal scenario 77 = 0. 



/i = 10 s GeV, 5=1, rj=0, taii/3=50, sgn(/()=+l, A =0 




200 400 600 800 1000 1200 1400 
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(c) Neutralino isomass contour in the plane 
(mo, m.1/2) for the non minimal scenario 
X] = 1. 



(d) Neutralino isomass contour in the plane 
(mo, m-ii 2) for the non minimal scenario 
77 = 2. 



Figure 7: Contour plots for iV = 2 Higgs models. 
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A /3-functions from the method of the effective po- 
tential 



In this appendix we collect the full set of renormalization group equations following the 
notation of jTJ] 

1 



Q j/i fabc 

dQ 



n d 
Q dQ^ ab 



d 



2\a 



Q dQ [m > b 



Q lQ i[l)a 



32vr 2 



Ag 2 [C(A)+C(B) + C(C)}f abc 



■(-X'o fa'bc + Xl fab'c + X° c f a bc' 

~ W[C(A) + C(B)]f, ab - (X a a '^ a , b + xU ah , 

1 



32tt 2 



32vr 2 



- X a c (m 2 \ - X c b (m 2 T c - 4f ade f bce (m 2 ) c d + 
-8£ (3 )a C( 4 C ) + 8g 2 [C(A)A 2 A + C(B)A 2 B ]5 a b 
-Ag 2 tr(T A m 2 )(T A ) a b 

+ Ag 2 [C(A)+C(B) + C(C)]^ )abc 

O / r fb'c'ds- I c ea'c'ds- , e ra'b'ds- 

— £ yjbcdj ?(3)afe'c' + JacdJ £(3)a'6c' + JabdJ £(3)a'b'. 

— Xa £(3)a'fcc — Xl C(3)a6'c ~ X° £(3)aftc' 

+4g 2 [C(A)A A + C(B)A B + C(C)A c }f abc 

1 r + Ag 2 [C(A) + C(B)]^ 2)ab - [X^ma'b + X^ (2)abf 



32tt 2 



2/abc/ C e ^(2)de — 4/ C e [ici)acd^eb + £(3)acd/^e&] 



+8( ? 2 [C(A)A A + C( J B)A i3 ]/i 



a6 



+ VC(A)e ( l )a " X^ (1)6 - 2f C Va fe £(2W 



327T 2 

-4^)^) - 4/ afec/ / d (m 2 )d 



(A.l) 
(A.2) 



(A.3) 



(A.4) 



(A.5) 



(A.6) 



The generalization of these equations to non simple gauge groups is very simple. Terms 
like g 2 [C(A) + C(B)] are to be replaced by a sum over the gauge groups (all simple 
components and U(l) factors). There is a gauge coupling constant and a group invariant 
C(R) for each group factor. 
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B Feynman diagrams analysis of the soft terms 



In this appendix we collect all the Feynman diagrams which renormalize the soft terms in 
(|2.1U|) . Next to each diagram we give its contribution to the coefficient of the respective 
/3-function. As one can easily check, in the absence of KK states, every diagram is loga- 
rithmically divergent in four dimensions. In presence of extra dimensions, those diagrams 
in which KK states can circulate contribute to the power law piece in ()2.32j) . This is the 
recipe we have adopted to obtain (l3~TUJ) - (|3~2UJ) and (j3~38|) - (j3~T7j) 

• Renormalization of the coefficient b in ()2.10|) 



t,b,r 



H 





(B.l) 



t,b,f 



H 



/ 



H 



2ii(3a t y t + 3a b y b + a T y T ) 



(B.2) 



\ 





(B.3) 



H 



A 





(B.4) 



H 
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Renormalization of the coefficient m 2 H , m 2 H in ()2.10|) 



t,b 



H f \ H no / \ 

-x- 6y 2 m 2 (B.5) 




t 

6a 2 t (B.6) 




A 

H H -§0?(A 11 ) 2 -6^(A 2a ) a (B.7) 

5 

H 

t 

- H — — ±* — — H - Qy 2 (m 2 + ml) (B.8) 



Jl-.-hA-.Jl -gftriYm 2 ) (B.9) 

5 



The last diagram is present only for the U(l) gauge group. 
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Renormalization of the coefficient m~ in (|2.10j) 




H 

f -*-- (2y 2 + 2y 2 )m 2 (B.fO) 

t 



H 

- — ( \-~ 2(a 2 + a 2 b ) (B.ff) 

t 



A 

/ ' -S,rAV'(/) (B.l-2) 

t 

t,H 

1 * 2y 2 (m 2 Hu + ml) + 2y 2 (m 2 Hd + mf) (B.13) 



t,H 

_? * -gltr(Ym 2 ) (B.14) 

5 




Similar diagrams hold for the b squark and the r slepton. 
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Renormalization of the coefficient a t , ay,, a T in ()2.f0|) 



\ t 



\ t 



t As / H 

-< at(3yt+Vb) 




\ 



t st 



t,b,r , _ 



v — ^ \ H 

t,b,r 



t,b,r 



S N 



t 



X S 



H — <f ' ' -W (y ^ A 33 + S^ 2 A 22 + |^ 2 A 11 



t — <f 1 ~Vt (^91 A 33 + 3s 2 A 22 + g, 2 A 11 

H ^ # 
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Similar diagrams hold for the other two trilinear couplings, a&, a T , the one for the squark 
b and the one for the slepton f. 
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